Defining and benchmarking open problems in single-cell analysis.

With the growing number of single-cell analysis tools, benchmarks are increasingly important to guide analysis and method development. However, a lack of standardisation and extensibility in current benchmarks limits their usability, longevity, and relevance to the community. We present Open Problems, a living, extensible, community-guided benchmarking platform including 10 current single-cell tasks that we envision will raise standards for the selection, evaluation, and development of methods in single-cell analysis.

Early life adversity shapes social subordination and cell type-specific transcriptomic patterning in the ventral hippocampus.

Adverse events in early life can modulate the response to additional stressors later in life and increase the risk of developing psychiatric disorders. The underlying molecular mechanisms responsible for these effects remain unclear. Here, we uncover that early life adversity (ELA) in mice leads to social subordination. Using single-cell RNA sequencing (scRNA-seq), we identified cell type-specific changes in the transcriptional state of glutamatergic and GABAergic neurons in the ventral hippocampus of ELA mice after exposure to acute social stress in adulthood. These findings were reflected by an alteration in excitatory and inhibitory synaptic transmission induced by ELA in response to acute social stress. Finally, enhancing the inhibitory network function through transient diazepam treatment during an early developmental sensitive period reversed the ELA-induced social subordination. Collectively, this study significantly advances our understanding of the molecular, physiological, and behavioral alterations induced by ELA, uncovering a previously unknown cell type-specific vulnerability to ELA.

Ex vivo tissue perturbations coupled to single-cell RNA-seq reveal multilineage cell circuit dynamics in human lung fibrogenesis.

Pulmonary fibrosis develops as a consequence of failed regeneration after injury. Analyzing mechanisms of regeneration and fibrogenesis directly in human tissue has been hampered by the lack of organotypic models and analytical techniques. In this work, we coupled ex vivo cytokine and drug perturbations of human precision-cut lung slices (hPCLS) with single-cell RNA sequencing and induced a multilineage circuit of fibrogenic cell states in hPCLS. We showed that these cell states were highly similar to the in vivo cell circuit in a multicohort lung cell atlas from patients with pulmonary fibrosis. Using micro-CT-staged patient tissues, we characterized the appearance and interaction of myofibroblasts, an ectopic endothelial cell state, and basaloid epithelial cells in the thickened alveolar septum of early-stage lung fibrosis. Induction of these states in the hPCLS model provided evidence that the basaloid cell state was derived from alveolar type 2 cells, whereas the ectopic endothelial cell state emerged from capillary cell plasticity. Cell-cell communication routes in patients were largely conserved in hPCLS, and antifibrotic drug treatments showed highly cell type-specific effects. Our work provides an experimental framework for perturbational single-cell genomics directly in human lung tissue that enables analysis of tissue homeostasis, regeneration, and pathology. We further demonstrate that hPCLS offer an avenue for scalable, high-resolution drug testing to accelerate antifibrotic drug development and translation.

Integrated multi-omics single cell atlas of the human retina.

Single-cell sequencing has revolutionized the scale and resolution of molecular profiling of tissues and organs. Here, we present an integrated multimodal reference atlas of the most accessible portion of the mammalian central nervous system, the retina. We compiled around 2.4 million cells from 55 donors, including 1.4 million unpublished data points, to create a comprehensive human retina cell atlas (HRCA) of transcriptome and chromatin accessibility, unveiling over 110 types. Engaging the retina community, we annotated each cluster, refined the Cell Ontology for the retina, identified distinct marker genes, and characterized cis-regulatory elements and gene regulatory networks (GRNs) for these cell types. Our analysis uncovered intriguing differences in transcriptome, chromatin, and GRNs across cell types. In addition, we modeled changes in gene expression and chromatin openness across gender and age. This integrated atlas also enabled the fine-mapping of GWAS and eQTL variants. Accessible through interactive browsers, this multimodal cross-donor and cross-lab HRCA, can facilitate a better understanding of retinal function and pathology.

A reversible state of hypometabolism in a human cellular model of sporadic Parkinson's disease.

Sporadic Parkinson's Disease (sPD) is a progressive neurodegenerative disorder caused by multiple genetic and environmental factors. Mitochondrial dysfunction is one contributing factor, but its role at different stages of disease progression is not fully understood. Here, we showed that neural precursor cells and dopaminergic neurons derived from induced pluripotent stem cells (hiPSCs) from sPD patients exhibited a hypometabolism. Further analysis based on transcriptomics, proteomics, and metabolomics identified the citric acid cycle, specifically the α-ketoglutarate dehydrogenase complex (OGDHC), as bottleneck in sPD metabolism. A follow-up study of the patients approximately 10 years after initial biopsy demonstrated a correlation between OGDHC activity in our cellular model and the disease progression. In addition, the alterations in cellular metabolism observed in our cellular model were restored by interfering with the enhanced SHH signal transduction in sPD. Thus, inhibiting overactive SHH signaling may have potential as neuroprotective therapy during early stages of sPD.

A Roadmap for a Consensus Human Skin Cell Atlas and Single-Cell Data Standardization.

Single-cell technologies have become essential to driving discovery in both basic and translational investigative dermatology. Despite the multitude of available datasets, a central reference atlas of normal human skin, which can serve as a reference resource for skin cell types, cell states, and their molecular signatures, is still lacking. For any such atlas to receive broad acceptance, participation by many investigators during atlas construction is an essential prerequisite. As part of the Human Cell Atlas project, we have assembled a Skin Biological Network to build a consensus Human Skin Cell Atlas and outline a roadmap toward that goal. We define the drivers of skin diversity to be considered when selecting sequencing datasets for the atlas and list practical hurdles during skin sampling that can result in data gaps and impede comprehensive representation and technical considerations for tissue processing and computational analysis, the accounting for which should minimize biases in cell type enrichments and exclusions and decrease batch effects. By outlining our goals for Atlas 1.0, we discuss how it will uncover new aspects of skin biology.

Synthetic enforcement of STING signaling in cancer cells appropriates the immune microenvironment for checkpoint inhibitor therapy.

Immune checkpoint inhibitors (ICIs) enhance anticancer immunity by releasing repressive signals into tumor microenvironments (TMEs). To be effective, ICIs require preexisting immunologically "hot" niches for tumor antigen presentation and lymphocyte recruitment. How the mutational landscape of cancer cells shapes these immunological niches remains poorly defined. We found in human and murine colorectal cancer (CRC) models that the superior antitumor immune response of mismatch repair (MMR)-deficient CRC required tumor cell-intrinsic activation of cGAS-STING signaling triggered by genomic instability. Subsequently, we synthetically enforced STING signaling in CRC cells with intact MMR signaling using constitutively active STING variants. Even in MMR-proficient CRC, genetically encoded gain-of-function STING was sufficient to induce cancer cell-intrinsic interferon signaling, local activation of antigen-presenting cells, recruitment of effector lymphocytes, and sensitization of previously "cold" TMEs to ICI therapy in vivo. Thus, our results introduce a rational strategy for modulating cancer cell-intrinsic programs via engineered STING enforcement to sensitize resistant tumors to ICI responsiveness.

Alveolar macrophages in early stage COPD show functional deviations with properties of impaired immune activation.

Despite its high prevalence, the cellular and molecular mechanisms of chronic obstructive pulmonary disease (COPD) are far from being understood. Here, we determine disease-related changes in cellular and molecular compositions within the alveolar space and peripheral blood of a cohort of COPD patients and controls. Myeloid cells were the largest cellular compartment in the alveolar space with invading monocytes and proliferating macrophages elevated in COPD. Modeling cell-to-cell communication, signaling pathway usage, and transcription factor binding predicts TGF-ß1 to be a major upstream regulator of transcriptional changes in alveolar macrophages of COPD patients. Functionally, macrophages in COPD showed reduced antigen presentation capacity, accumulation of cholesteryl ester, reduced cellular chemotaxis, and mitochondrial dysfunction, reminiscent of impaired immune activation.

Primary cilia and SHH signaling impairments in human and mouse models of Parkinson's disease.

Parkinson's disease (PD) as a progressive neurodegenerative disorder arises from multiple genetic and environmental factors. However, underlying pathological mechanisms remain poorly understood. Using multiplexed single-cell transcriptomics, we analyze human neural precursor cells (hNPCs) from sporadic PD (sPD) patients. Alterations in gene expression appear in pathways related to primary cilia (PC). Accordingly, in these hiPSC-derived hNPCs and neurons, we observe a shortening of PC. Additionally, we detect a shortening of PC in PINK1-deficient human cellular and mouse models of familial PD. Furthermore, in sPD models, the shortening of PC is accompanied by increased Sonic Hedgehog (SHH) signal transduction. Inhibition of this pathway rescues the alterations in PC morphology and mitochondrial dysfunction. Thus, increased SHH activity due to ciliary dysfunction may be required for the development of pathoetiological phenotypes observed in sPD like mitochondrial dysfunction. Inhibiting overactive SHH signaling may be a potential neuroprotective therapy for sPD.

Ly6D+Siglec-H+ precursors contribute to conventional dendritic cells via a Zbtb46+Ly6D+ intermediary stage.

Plasmacytoid and conventional dendritic cells (pDC and cDC) are generated from progenitor cells in the bone marrow and commitment to pDCs or cDC subtypes may occur in earlier and later progenitor stages. Cells within the CD11c+MHCII-/loSiglec-H+CCR9lo DC precursor fraction of the mouse bone marrow generate both pDCs and cDCs. Here we investigate the heterogeneity and commitment of subsets in this compartment by single-cell transcriptomics and high-dimensional flow cytometry combined with cell fate analysis: Within the CD11c+MHCII-/loSiglec-H+CCR9lo DC precursor pool cells expressing high levels of Ly6D and lacking expression of transcription factor Zbtb46 contain CCR9loB220hi immediate pDC precursors and CCR9loB220lo (lo-lo) cells which still generate pDCs and cDCs in vitro and in vivo under steady state conditions. cDC-primed cells within the Ly6DhiZbtb46- lo-lo precursors rapidly upregulate Zbtb46 and pass through a Zbtb46+Ly6D+ intermediate stage before acquiring cDC phenotype after cell division. Type I IFN stimulation limits cDC and promotes pDC output from this precursor fraction by arresting cDC-primed cells in the Zbtb46+Ly6D+ stage preventing their expansion and differentiation into cDCs. Modulation of pDC versus cDC output from precursors by external factors may allow for adaptation of DC subset composition at later differentiation stages.

The discovAIR project: a roadmap towards the Human Lung Cell Atlas.

The Human Cell Atlas (HCA) consortium aims to establish an atlas of all organs in the healthy human body at single-cell resolution to increase our understanding of basic biological processes that govern development, physiology and anatomy, and to accelerate diagnosis and treatment of disease. The Lung Biological Network of the HCA aims to generate the Human Lung Cell Atlas as a reference for the cellular repertoire, molecular cell states and phenotypes, and cell-cell interactions that characterise normal lung homeostasis in healthy lung tissue. Such a reference atlas of the healthy human lung will facilitate mapping the changes in the cellular landscape in disease. The discovAIR project is one of six pilot actions for the HCA funded by the European Commission in the context of the H2020 framework programme. discovAIR aims to establish the first draft of an integrated Human Lung Cell Atlas, combining single-cell transcriptional and epigenetic profiling with spatially resolving techniques on matched tissue samples, as well as including a number of chronic and infectious diseases of the lung. The integrated Human Lung Cell Atlas will be available as a resource for the wider respiratory community, including basic and translational scientists, clinical medicine, and the private sector, as well as for patients with lung disease and the interested lay public. We anticipate that the Human Lung Cell Atlas will be the founding stone for a more detailed understanding of the pathogenesis of lung diseases, guiding the design of novel diagnostics and preventive or curative interventions.

Mapping single-cell data to reference atlases by transfer learning.

Large single-cell atlases are now routinely generated to serve as references for analysis of smaller-scale studies. Yet learning from reference data is complicated by batch effects between datasets, limited availability of computational resources and sharing restrictions on raw data. Here we introduce a deep learning strategy for mapping query datasets on top of a reference called single-cell architectural surgery (scArches). scArches uses transfer learning and parameter optimization to enable efficient, decentralized, iterative reference building and contextualization of new datasets with existing references without sharing raw data. Using examples from mouse brain, pancreas, immune and whole-organism atlases, we show that scArches preserves biological state information while removing batch effects, despite using four orders of magnitude fewer parameters than de novo integration. scArches generalizes to multimodal reference mapping, allowing imputation of missing modalities. Finally, scArches retains coronavirus disease 2019 (COVID-19) disease variation when mapping to a healthy reference, enabling the discovery of disease-specific cell states. scArches will facilitate collaborative projects by enabling iterative construction, updating, sharing and efficient use of reference atlases.

Benchmarking atlas-level data integration in single-cell genomics.

Single-cell atlases often include samples that span locations, laboratories and conditions, leading to complex, nested batch effects in data. Thus, joint analysis of atlas datasets requires reliable data integration. To guide integration method choice, we benchmarked 68 method and preprocessing combinations on 85 batches of gene expression, chromatin accessibility and simulation data from 23 publications, altogether representing >1.2 million cells distributed in 13 atlas-level integration tasks. We evaluated methods according to scalability, usability and their ability to remove batch effects while retaining biological variation using 14 evaluation metrics. We show that highly variable gene selection improves the performance of data integration methods, whereas scaling pushes methods to prioritize batch removal over conservation of biological variation. Overall, scANVI, Scanorama, scVI and scGen perform well, particularly on complex integration tasks, while single-cell ATAC-sequencing integration performance is strongly affected by choice of feature space. Our freely available Python module and benchmarking pipeline can identify optimal data integration methods for new data, benchmark new methods and improve method development.

Diet-induced alteration of intestinal stem cell function underlies obesity and prediabetes in mice.

Excess nutrient uptake and altered hormone secretion in the gut contribute to a systemic energy imbalance, which causes obesity and an increased risk of type 2 diabetes and colorectal cancer. This functional maladaptation is thought to emerge at the level of the intestinal stem cells (ISCs). However, it is not clear how an obesogenic diet affects ISC identity and fate. Here we show that an obesogenic diet induces ISC and progenitor hyperproliferation, enhances ISC differentiation and cell turnover and changes the regional identities of ISCs and enterocytes in mice. Single-cell resolution of the enteroendocrine lineage reveals an increase in progenitors and peptidergic enteroendocrine cell types and a decrease in serotonergic enteroendocrine cell types. Mechanistically, we link increased fatty acid synthesis, Ppar signaling and the Insr-Igf1r-Akt pathway to mucosal changes. This study describes molecular mechanisms of diet-induced intestinal maladaptation that promote obesity and therefore underlie the pathogenesis of the metabolic syndrome and associated complications.

Single-cell meta-analysis of SARS-CoV-2 entry genes across tissues and demographics.

Angiotensin-converting enzyme 2 (ACE2) and accessory proteases (TMPRSS2 and CTSL) are needed for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) cellular entry, and their expression may shed light on viral tropism and impact across the body. We assessed the cell-type-specific expression of ACE2, TMPRSS2 and CTSL across 107 single-cell RNA-sequencing studies from different tissues. ACE2, TMPRSS2 and CTSL are coexpressed in specific subsets of respiratory epithelial cells in the nasal passages, airways and alveoli, and in cells from other organs associated with coronavirus disease 2019 (COVID-19) transmission or pathology. We performed a meta-analysis of 31 lung single-cell RNA-sequencing studies with 1,320,896 cells from 377 nasal, airway and lung parenchyma samples from 228 individuals. This revealed cell-type-specific associations of age, sex and smoking with expression levels of ACE2, TMPRSS2 and CTSL. Expression of entry factors increased with age and in males, including in airway secretory cells and alveolar type 2 cells. Expression programs shared by ACE2+TMPRSS2+ cells in nasal, lung and gut tissues included genes that may mediate viral entry, key immune functions and epithelial-macrophage cross-talk, such as genes involved in the interleukin-6, interleukin-1, tumor necrosis factor and complement pathways. Cell-type-specific expression patterns may contribute to the pathogenesis of COVID-19, and our work highlights putative molecular pathways for therapeutic intervention.

Current best practices in single-cell RNA-seq analysis: a tutorial.

Single-cell RNA-seq has enabled gene expression to be studied at an unprecedented resolution. The promise of this technology is attracting a growing user base for single-cell analysis methods. As more analysis tools are becoming available, it is becoming increasingly difficult to navigate this landscape and produce an up-to-date workflow to analyse one's data. Here, we detail the steps of a typical single-cell RNA-seq analysis, including pre-processing (quality control, normalization, data correction, feature selection, and dimensionality reduction) and cell- and gene-level downstream analysis. We formulate current best-practice recommendations for these steps based on independent comparison studies. We have integrated these best-practice recommendations into a workflow, which we apply to a public dataset to further illustrate how these steps work in practice. Our documented case study can be found at https://www.github.com/theislab/single-cell-tutorial This review will serve as a workflow tutorial for new entrants into the field, and help established users update their analysis pipelines.